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Abstract 

Future power grids will offer enhanced controllability due to the increased avail¬ 
ability of power flow control units (FACTS). As the installation of control units in the 
grid is an expensive investment, we are interested in using few controllers to achieve 
high controllability. In particular, two questions arise: How many flow control buses 
are necessary to obtain globally optimal power flows? And if fewer flow control buses 
are available, what can we achieve with them? 

Using steady state IEEE benchmark data sets, we explore experimentally that al¬ 
ready a small number of controllers placed at certain grid buses suffices to achieve 
globally optimal power flows. We present a graph-theoretic explanation for this behav¬ 
ior. To answer the second question we perform a set of experiments that explore the 
existence and costs of feasible power flow solutions at increased loads with respect to 
the number of flow control buses in the grid. We observe that adding a small number of 
flow control buses reduces the flow costs and extends the existence of feasible solutions 
at increased load. 


1 Introduction 

The central task of any electrical power infrastructure is the reliable and cost-efficient supply 
of electrical energy to industry and population on a national or even continental scale. Future 
power grids and their usage are subject to fundamental changes due to the shift towards 
renewable distributed energy production and the installation of new power flow control 
units, which offer increased control, but make the grid operation more demanding. Not 
only do these changes lead to a much larger number of independent power producers (IPP), 
which are highly distributed in the network, but they also cause very different patterns of 
energy flow. For example, regions with off-shore wind farms may sometimes produce enough 
energy to supply remote consumers, but at other times they are consumers themselves. In 
particular, this may require long-distance energy transmission and frequent flow direction 
changes. Most of the existing power grids, however, were not designed for such transmission 
patterns. The current strategy to cope with these changes is to either extend the grid 
with additional transmission lines, or to install advanced control units to facilitate better 
utilization of the existing infrastructure. 
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In this paper, we consider the latter option and study the advantageous effects of making 
selected buses of a power grid to controllable, both in terms of the minimum number of 
controllable buses needed for achieving maximum flow control and in terms of the operation 
costs and the existence of feasible power flows at critical line capacities. 

In abstract terms, we assume that a flow-control bus is able to flexibly distribute the 
entire power flow at this bus among the incident edges, as long as Kirchhoff’s current law 
(or the flow-conservation property) is satisfied, i.e., the in-flow to the bus equals its out¬ 
flow. These flow control buses can be realized using power electronics devices known as 
flexible AC transmission systems (FACTS), which are a class of power systems that have 
the capabilities to control various electrical bus parameters II T^. More specifically, since 


we are interested in controlling the real power flow on the branches incident to a particular 
bus, we can realize our flow control buses by installing on each (but one) incident branch a 
unified power flow controller (UPFC), which is a FACTS that is able to control the voltage 
magnitude and angle and consequently has control of the real and reactive power flow on 
the particular branch 
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One of the most important tasks in operating a power grid is to decide the energy pro¬ 
duction of each power generator such that supply and consumption are balanced and the 
resulting power flow does not exceed the thermal limits of the power lines. Among all solu¬ 
tions we are interested in one that minimizes the total energy production and transmission 
cost. This is called Economic Dispatch Problem (EDP). The main approach for solving 
this problem in power grids without FACTS is the optimal power flow (OPF) method, a 
numerical method that was introduced in 1962 by Carpentier and has subsequently been 
refined and generalized, see the recent surveys by Frank et al. [^[^. However, OPF is not 
designed for hybrid power grids with flow control buses and cannot exploit the extended 
flow control possibilities to obtain globally optimal solutions. 

Hence, we propose in Section]^ a new hybrid DC-based model for power flows in power 
grids that combine traditional grid buses with some flow-control buses. In order to answer 
our questions on the effects of installing flow control buses, we solve the EDP in our hybrid 
model using a linear programming (LP) formulation. Our LP combines a standard graph- 
theoretical network flow model, which already includes Kirchhoff’s current law at all buses, 
with additional constraints for Kirchhoff’s voltage law in those parts of the grid that are 
not equipped with flow control buses. Thus we are able to obtain electrically feasible power 
flows that minimize, similarly to OPF, the overall flow costs in terms of generation and 
transmission costs. 

Using the well-known IEEE power systems test cases, we performed simulation exper¬ 
iments related to two key questions, which take into account that the FACTS needed for 
realizing our flow control buses in reality constitute a significant and expensive investment 
and hence their number should be as small as possible. 

1. How many flow control buses are necessary to obtain globally optimal power flows and 
which buses need to be controlled? 

2. If the number of available flow control buses is given, do we still see a positive effect on 
the flow costs and on the operability of the grid when approaching its capacity limits? 

In Sectionj^we address the first question. In our experiments we determine the minimum 
number of flow control buses necessary to achieve the same solution quality as in a power 
grid in which each bus is controllable and which clearly admits an upper bound on what 
can be achieved with the network topology. Interestingly, it turns out that a relatively 
small number of flow control buses are sufficient for this. In fact, we can prove a theorem 
stating a structural graph-theoretic property, which, if met by the placement of flow control 
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buses, implies the optimality of the power flow and serves as a theoretical explanation of the 
observed behavior. Sectiondeals with the second question of operating a power grid close 
to its capacity limits, which becomes increasingly relevant as the consumption of electrical 
energy grows faster than the grid capacities. Our experiments indicate that installing few 
flow control buses in a power grid is sufficient not only to achieve lower costs compared to 
an OPF solution, but also allows to operate the grid at capacities for which no feasible OPF 
solution exists any more. 


2 Related Work 


With the increasing availability and technological advancement of FACTS researchers began 
to study the possible benefits of their installation in power grids from different perspectives. 

From an economic perspective, it is of interest to support investment decisions in power 
grid expansion planning by considering alternative investment strategies that either focus 
on new transmission lines or allow mixed approaches including FACTS placement. Blanco 
et al. present a least-squares Monte-Carlo method for evaluating investment strategies 
and argue that FACTS allow for a more flexible, mixed strategy that fares better under 
uncertainty. Tee and flic present an optimal decision-making framework for comparing 
investment decisions, including FACTS 
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From the perspective of operating a power grid, the main question is how many and 
where FACTS should be placed in order to optimize a certain criterion. Cai et al. 
propose and experimentally evaluate a genetic algorithm for allocating different types of 
FACTS in a power grid in order to optimally support a deregulated energy market. Gerbex 
et al. and Ongsakul and Jirapong 17 study the placement of FACTS with the goal 
of increasing the amount of energy that can be transferred. Gerbex et al. present a 
genetic algorithm that optimizes simultaneously the energy generation costs, transmission 
losses, line overload, and the acquisition costs for FAGTS. Ongsakul and Jirapong 17 


use 


evolutionary programming to place FACTS such that the total amount of energy that can 
be transferred from producers to consumers is maximized; in contrast to our setting, they 
may also increase the demands of consumers arbitrarily. In contrast to these heuristic 
approaches Lima et al. 15 use mixed-integer linear programming to optimally increase the 


loadability of a system by placing FACTS subject to limits on their number or cost. Similar 
to our approach, they do not distinguish different types of FACTS but rather assume “ideal” 
FACTS that can control all transmission parameters of a branch. In contrast to our work, 
they focus only on loadability and do not consider generation costs and line losses. 

All related work mentioned so far considers the DC model for electrical networks as 
an approximation to the AC model and aims at providing a preliminary step in an actual 
planning process, where this approximation is sufficient. There are also a few attempts to 
solve the placement problem for FACTS in the more realistic but also more complicated AC 
model. Sharma et al. 20 develop an evaluation whether transmission lines are critical and 


propose to place FACTS at critical lines in order to improve voltage stability in the grid. 


Ippolito and Siano 14 present a genetic algorithm for FACTS placement in AC networks and 


experimentally evaluate it in a case study. In contrast to these heuristic approaches, Farivar 
and Low observe exact OPF evaluation in a relaxed AC-model. In this context, they 
place phase shifters to exploit structural characteristics that are similar to our approach. 
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3 Preliminaries 


In this section we recall some basic notions from graph theory. Although, for technical 
reasons, the graphs we use for modeling power grids are directed, when considering the 
topology of the network, we always consider the underlying undirected graph. Thus, in the 
following let G be an undirected graph. 

The graph G is connected if it contains a path between any two vertices. A connected 
component of G is a maximal connected subgraph of G (maximal with respect to inclusion). 
A cactus is a graph where every edge is contained in at most one cycle. A forest is a graph 
that does not contain a cycle. 

A cutvertex is a vertex of a graph whose removal increases the number of connected com¬ 
ponents. A hiconnected component is a maximal subgraph that does not have a cutvertex. 
Note that a biconnected component of a forest is either trivial in the sense that it consists 
of a single vertex, or it consists of a single edge. Similarly, a biconnected component of a 
cactus is trivial, a single edge, or a cycle. 

A feedback set of G = {V, E) with respect to a class of graphs ^ is a set of vertices F C V 
such that G — F G G- We will only be interested in feedback sets with respect to forests and 
cacti. The former is also called feedback vertex set. Naturally, one is interested in finding a 
set F that is as small as possible. 


4 Model 


In this section we introduce three graph-theoretic flow models for optimal power flows. Our 
models are based on the DC power grid model 12,21 25 , which is commonly used as an 


approximation of AC grids [l8p^ . We model a power grid as a graph G = {V, E), where V is 
the set of vertices and E C (^ is the set of edges. The vertices represent the buses, some of 
which may be special generator and consumer buses, and the edges represent the branches, 
which may be transmission lines between the incident buses or transformers. There is a 
subset Vg GV oi the vertices that represents generator buses. Each generator g GVq has 
a maximum supply Xg G and is equipped with a convex cost function 7 g > 0 that is 
assumed to be piecewise linear with 


')g[x) = max{aja; -I- Ci \ (oi, cf) G Fg), (1) 

where Fg is the set of all piecewise linear functions of 7 g and < a^+i. Further, there is a 
subset Vc GV\Vg consumer buses. Each consumer u G Vc has a power demand G K. 

Each branch e G E has a thermal limit, which is modeled as a capacity function c: A —>■ M 
restricting the real power. Further, each branch causes a certain loss of power depending 
on the physical branch parameters and the actual power flow on the branch. These losses 
are again approximated as a convex, piecewise linear function for each edge e G E with 


£e{x) = max{aia; -I- c* | {at, cf) G Fe}, 


( 2 ) 


where F^ is the set of all piecewise linear functions of and Oi < a^+i. 

A flow f in the power grid G is a function /: E x E —)■ K with the property that 


f{u, v) = -f{v, u) yu,v GV 


( 3 ) 
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For every vertex u in G, we define its net out-flow 

fnet{u)= f{u,v). 

{u^v}gE 

For a flow /, we further define two types of costs, the generator costs 

=»(/) = E 7g(/net(5)) 
geVa 


and the line losses 

Gif)= Y ■ 

{u^v}£E 

To obtain the overall cost for the flow /, we weight these two terms as 

cxif) = X ■ Cgif) + (1 - A) • cflf) (4) 

where A G [0,1]. Our goal is to minimize this objective function in several different power 
flow models. 


4.1 Power Flow Models 

The most basic model is the flow model, where / has to satisfy the following constraints. 


-c(e) < f{u,v) < c(e) 

fnet{v) = 0 

/net('c) = dy 

0 E fnetii^) E 


Ve = {u, v} £ E 

(5) 

u e E \ (Eg U Eg) 

(6) 

w G Eg 

(7) 

u G Eg 

(8) 


We call a flow satisfying these constraints feasible. Equation ([^ models the thermal limits or 
real power capacities of all branches and is called capacity constraint. Equation ®> models 
that vertices that are neither generators nor consumers have zero net out-flow and is called 
flow conservation constraint. Equation Q models that all consumer demands are satisfied 
and is called consumer constraint. Finally, Equation <§ models that all generators respect 
their production limits and is called generator constraint. 

The flow model neglects some physical properties of electrical flows, in particular Kirch- 
hoff’s voltage law. Thus, the computed power flows can only be applied to power grids 
where every vertex is a control vertex. In contrast, the electrical flow model, e.g., according 
to Zimmerman et al. 25 , models the power flow via the same set of constraints as the 


flow model, but additionally requires the existence of a suitable voltage angle assignment 
0: E —)■ K such that for each branch {u,z;} the following equation holds 


f{u,v) = B{u,v){e{u) - e{v)). 


(9) 


Here B{u,v) is the susceptance of the branch {u,v). This is equivalent to restricting the 
model to feasible flows that also satisfy Kirchhoff’s voltage law, or, in other words, no flow 
control buses are used. This yields a model that matches the situation in the traditional 
power grids existing today. We call a feasible flow / electrically feasible if there exists a 
voltage angle assignment 0 satisfying 
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case 

nb 

nl 

ng 

Pd 

case6 

6 

11 

3 

210.00 

case9 

9 

9 

3 

315.00 

casel4 

14 

20 

5 

259.00 

caseSO 

30 

41 

6 

189.20 

case39 

39 

46 

10 

6254.23 

case57 

57 

78 

7 

1250.80 

casellS 

118 

179 

54 

4242.00 


Table 1: IEEE benchmark set with n&, nl, ng and pd representing number of buses, number 
of transmission lines, number of generators and total power demand, respectively. 


Recall from the introduction that flow control buses can be technically realized by 
UPFCs, which is a FACTS. Ideal FACTS as introduced by Griffin et al. 10 are often 


used to simplify the modeling of FACTS by using a linear model and assuming a complete 
and independent control of the real and reactive power. Our flow control buses are ideal 
FACTS that control the power flow to all incident edges. The flow model—in contrast to the 
electrical model—assumes flow control buses at each vertex, whereas the electrical model 
assumes no immediate control of the power flow. Instead, the grid is balanced by changing 
the generator outputs only. In the following we propose a hybrid model that combines the 
flow model and the electrical flow model in order to handle power grids with flow control 
buses at a subset of selected vertices. 

Let F C V he a subset of vertices of G. We denote hy Gp the power network obtained 
from G by considering all vertices in F as flow control buses. We call any subgraph G" = 
G\V'] induced by a subset V' C V \ F oi the vertices without controllers a native power 
grid of G. A flow of G is electrically feasible for a native power grid G' C G if there exists 
a voltage angle assignment 0: if —>• K such that every edge in G' satisfies Equation . In 
this case we call 0 feasible (voltage) angle assignment for G'. 

A feasible flow / is electrically feasible for Gp if and only if / is electrically feasible for 
the maximal native power grid G — F = G\y \ F]. Intuitively, this models the fact that a 
power flow in Gf must be a feasible flow and that it satisfies the second Kirchhoff law in 
the maximum native power grid. 

Obviously, if F C F' and / is an electrically feasible flow for Gf, then / is also electrically 
feasible for Gpi. Hence the minimum value of the cost Cx does not increase when adding 
more flow control buses. 

We note that each of the models can easily be expressed as a linear program (LP), and 
thus in all three models an optimal solution can be computed efficiently [^; see Appendix [A| 
However, the flow model can be reduced to a special minimum cost network flow problem, 
for which efficient exact optimization algorithms exist . We describe this reduction in the 
following. 


4.2 Reduction to MinCostFlow 

Let N = {G = (P, F); (s, t); c; a) be an s-t flow network consisting of a directed (multi-) 
graph G, two dedicated source and sink vertices s,t &V, edge capacities c : E ^ Kg , and 
edge costs a : E ^ K)}". A flow / in A^ is a function f : E ^ K)}" and it is called feasible if 
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( 10 ) 


it satisfies the capacity constraint ([^ and a flow conservation constraint similar to ([^: 

f{u,v)- Y /(^^>'«)=0 V-y e {s,t} 

{u,v)£E {v,u)£E 

The value \f\ of a flow / is the total flow from s to t, i.e., |/| = E{ut)eE f{u,t) = 
E(s u)g-E^ feasible flow / with maximum value is called a maximum flow in N. 
For a given flow value b the min-cost s-t flow problem is to find a feasible flow / of value 
I/I = b such that the cost cat)/) = ’ /(®) minimized. 

In order to transform the graph G = (V, E) of a power grid into an s-t flow network N, 
we first add a new source vertex s and a new sink vertex t to V. Each generator g G Vg 
is connected by a directed edge {s,g) with capacity c{s,g) = Xg to the source s. Each 
consumer u € Vc is connected by a directed edge (u, t) with capacity c(it, t) = to the 
sink t. Further, we replace each original undirected edge {u,y} € E by two directed copies 
(u, v) and {v, u), whose capacities c(m, v) = c(v, u) are given by their common thermal limit 
c({u,y}). 

Next, we define the edge costs. It is well known that a convex, piecewise linear edge cost 
function h can easily be modeled in a flow network by replacing the respective edge [u, v) 
with as many copies as the linear pieces of the cost function. The edge capacities are defined 
by the differences between consecutive breakpoints of h and sum up to c(u, y); the individual 
costs correspond to the costs as defined by the linear pieces between the breakpoints. Thus 
for ease of presentation we refrain from explicitly modeling convex piecewise linear cost 
functions in N. We rather assume that the flow cost is given for each edge (u, y) {u ^ s, 
V ^ t) hy the weighted loss function z\{{u,v), f ) = (1 — A) • •^{«,«}(/(«, y)), where A is the 
weight parameter of Equation Q. The edges {s,g) from the source s to a generator g have 
cost z\{{s,g), f) = X ■ yg{f{s,g)) and the edges incident to the sink t have cost 0. Then 
the objective function to be minimized is z\{f) = Ee^E -^>(6, /). Finally, we set the target 
flow value b to the total demand EueVc consumers. By construction, every feasible 

minimum-cost flow in TV is a feasible minimum-cost flow in the underlying power grid G 
and vice versa. 


5 Placing Flow Control Buses 

In this section we seek to answer the question how many flow control buses are necessary to 
obtain a globally optimal solution. Recall that the flow model is a relaxation of the physical 
model and uses fewer constraints. Therefore, optimal solutions in the flow model are at 
least as good as in the physical model. 

Given a power grid G = (V, E), we say that making the vertices in F flow control buses 
achieves full control if the objective value of an optimal energy flow for the grid Gp is the 
same as the objective value of an optimal solution in the flow model (or equivalently in the 
hybrid network Gy, where every vertex is a flow control bus). Our experiments indicate 
that in the IEEE instances often a small fraction of the vertices is sufficient to achieve full 
control. Afterwards we give a graph-theoretical explanation of this behavior. 
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(a) Total costs and losses for the IEEE bench- (b) Relative costs and losses for the IEEE bench¬ 
mark data with 30 buses, where the square cross mark data sets, 
marks the solution computed by OFF. 

Figure 1: Trade-off of generator costs and losses normalized to the maximum generator cost 
(A = 0) and the maximum loss (A = 1) as A varies from 0 to 1. 


5.1 Experiments 

For our evaluation we use the IEEE benchmark data set^ shown in Table [T] There each 
case is named according to the number of buses nb. The number of generators and the 
number of edges are denoted ng and nl, respectively. 

To obtain piece-wise linear functions for generator costs and line losses, we simply sample 
the cost functions using a specified number of sampling points. Note further that our 
approach requires convex cost functions, but this is fine in practice 23 ; in particular the 
functions are convex for the IEEE benchmark instances. 

We performed our experiments on an AMD Opteron processor 6172 running openSUSE 
12.2. Our implementation is written in Python 2.7.3 and uses PYP0WEB[^ a Python 
port of MATPOWER 24 for computing OPE solutions. For computing solutions and 
minimizing the number of control buses in our hybrid model we use the (integer) linear 
programming solver Gurobi 6.0. 10 

First, we observe that the value of A, which controls the weighting of costs and losses 
in the objective value has a significant effect on the objective values of generator costs and 
line losses. Figure la shows the trade-off for the IEEE instance case30 (the plots for the 
other instances can be found in Appendix H- The OPE solution, which ignores losses, is 
typically at the far end of the spectrum with high losses and is comparable to our solution 
with A = 1. As can be seen in Figure lb where the costs and losses are normalized to the 
maximum cost and the maximum loss per instance, the same trade-off behavior is present 
in all instances. It thus makes sense to allow the operator of a power grid to choose the 
value of A in order to model the true operation costs. 

On the other hand, it may then be the case that the number of flow control buses to 


^data sources http://www.pserc.cornell.edu/matpower/ and http://www.ee.washington.edu/ 
research/pstca/ 

https : //pypi .python. org/pypi/PYPOWER/4.0.0 
^ WWW.gurobi.com 



























Figure 2: Relative number of controllers for achieving full control in the IEEE instances as 
A varies from 0 to 1. 


achieve full control of the network varies depending on the choice of A. Figure shows for 
different values of A the relative number of control vertices necessary to achieve full control 
in each of the instances. In most cases less than 15% of all buses need to be controllers 
to achieve full control. For the cases with 6 buses and 14 buses this percentage is slightly 
bigger, which is mainly an artifact stemming from the small total size. As can be seen, the 
required number of units is relatively stable but drops to zero for A = 1, i.e., when only the 
generator costs are considered. This is due to the fact that all IEEE instances have basically 
unlimited line capacities and thus do not restrict the possible flows. 

In order to make a useful prediction on the number of vertices required for full control 
that applies to all choices of A, in the following we take for each instance the maximum of 
the smallest possible number of vertices to achieve full control over all values of A and refer 
to this as the number of vertices for achieving full control of the instance. This conservative 
choice ensures that the numbers we compute are certainly an upper bound for achieving full 
control, independent of the actual choice of A. 

5.2 Structure of Optimal Solutions 

As we have seen in our experimental evaluation, often a small number of flow control buses is 
sufficient to ensure that solutions in the hybrid model are the same as in the flow model. In 
the following we provide a theoretical explanation of this property and link it to structural 
properties of power grids. Earivar and Low give similar structural results on spanning 
trees, but using a different model. 

A first observation is that flow control buses influence all incident edges. Thus, if every 
edge is incident to a flow control bus, i.e., the set F is a vertex cover of G, no edge in 
the network is affected by constraint Then the flow model and the hybrid model are 
equivalent and full control is achieved. However, it is generally not true that power grids 
admit small vertex covers; as shown in Figure all instances require more than 40% of 
their vertices for a vertex cover. In the following we show a much stronger result, namely 
that it suffices for becoming independent of Equation § that the native power grid G — F 
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Figure 3: Comparison of the number of vertices which need to be removed from the network 
to get a Tree (Feedback Vertex Set) or a cactus, with the worst number of controller to have 
full control in the network. 


is an acyclic network. Moreoever, if A = 1, (line losses are neglected) and edge capacities 
are ignored, it even suffices that G — F is a. so-called cactus graph, in which every edge is 
part of at most one cycle. 

Lemma 5.1. Let H = {V,E) he a native power grid and let v be a vertex whose removal 
disconnects H into connected components with vertex sets Ci,... ,Ck- Then a flow f is 
electrically feasible for H if and only if it is electrically feasible for F[i = H[Ci U {u}] for 
i = 1,... ,k. 

Proof. Clearly, if 0 is a feasible voltage angle assignment for H, then its restriction to 
CiU{v} is a feasible angle assignment for Hi. Conversely, assume that 0.^ is a feasible angle 
assignment for Hi. Dehne 0' = 0i — Qi{v). Since for every edge in Hi the voltage angles 
of the endpoints are changed by the same value, 0' is a feasible voltage angle assignment 
for Hi. Further, 0((u) = 0 for every Hi, which means that the function 0: V —)■ K, where 
0(u) !-)• 0((it) for u G Ci is well-defined. Note that the restriction of 0 to any of the Hi 
coincides with 0'. Since every edge of H belongs to exactly one of the Hi, it follows that 0 
is a feasible voltage angle assignment for H. □ 

Iteratively applying Lemma |5 .1 1 yields the following. 

Corollary 5.2. A flow in a native power grid is electrically feasible if and only if it is 
electrically feasible for each hiconnected component of the power grid. 

We observe that if G — F is a forest, then each biconnected component H consists of a 
single edge {u,v}. Then 0(u) = f{u,v)/B{u,v) and 0(u) = 0 are feasible voltage angles 
for any flow f in B. Thus we conclude with the following 

Theorem 5.3. Let H be a native power grid that is a forest. Then every flow f is electrically 
feasible on H. 
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Thus, when F is a feedback vertex set of G, i.e., G — F is a forest, then every flow on 
G is electrically feasible for G — F, and thus any feasible flow for Gp is electrically feasible 
for Gp- It follows that the flow model and the hybrid model are equivalent in this case. In 
particular, whenever F is a feedback vertex set, instead of solving the LP for the hybrid 
model, we can rather assume the flow model and compute an optimal solution using a, 
potentially more efficient, flow algorithm. It follows from Theorem |5.3| that this solution is 
optimal also in the hybrid model. 

Figure shows for each of our instances the relative number of vertices necessary to 
obtain a vertex cover, a feedback vertex set with respect to forests, and the number of 
vertices necessary to obtain full control. In all instances a vertex cover is two to three times 
larger than a feedback vertex set (for forests) and the vertex set necessary for full control. 
Comparing the relative number of controllers for full control with the size of a feedback 
vertex sets shows that the number to get an optimal placement is in many cases smaller 
than the size of a feedback vertex set. Thus, in the optimal solutions, the native power grid 
does not always represent a forest, but can also include cycles. A closer inspection showed 
that this is in particular the case for instances that are operated far from their capacity 
limits. 

We now consider what happens when cycles exist in a native power grid. To this end, 
we start with the simplest case of a power grid that consists of a single cycle G. We say 
that two flows / and f' on a network G = (V, E) are equivalent if for each vertex v G V we 
have fnetiv) = f^iv). 

Lemma 5.4. Let G be a native power grid that is a cycle. For every flow f there exists a 
unique equivalent flow f that is electrically feasible for G. 

Proof. Let z)i,... ,u„ be the vertices of G as they occur along the cycle, i.e., f{vi,Vj) = 0 
unless Vi and Vj are neighbors on the cycle. Assume we wish to change the amount of 
flow from vi to V 2 by a fixed amount A and obtain an equivalent flow. The net out-flow 
conservation at the vertices then uniquely determines the change of flow along the remaining 
edges. Hence, every flow /' equivalent to / is obtained from / by choosing some amount A 
and setting f{vi,Vi+i) = f{vi,Vi+i) + A and f{vi+i,Vi) = f{vi+i,Vi)-A, where Vn+i = vi. 

Now the existence of a suitable offset A and the associated feasible voltage angles can 
be expressed as a linear system of equations. Namely, for edge {vi, Ui+i), f = 1,..., n (again 
using Vn+i = vi), we have the equation 


B{vi,Vi+i) ■ Q{vi) - B{vi, Vi+i) ■ ©(r^i+i) - A = f{vi, Vi+i). 


It is readily seen that the n equations are linearly independent, and hence a solution ex¬ 
ists. Moreover, dividing each of the equations by B(vi,Vipi) and summing them up yields 
— ^"^ 1 1/F(ui, Ui+i)A = X]r=i^*- 1 - 1)1 which shows that the value A is 
uniquely determined. □ 


Note however, that the equivalent flow /' whose existence is guaranteed by Lemma 5.4 
does not necessarily satisfy the capacity constraints. Also the evaluation of f in terms 
of line losses may change. If neither of these is a limiting factor, e.g., if A = 1 and line 
capacities are sufficiently large, we can show a stronger version of Theorem |5.3[ Recall that 
a cactus is a graph where every edge belongs to at most one cycle. 


Theorem 5.5. Let Gp be a power grid with flow control buses at the vertices in F such 
that the maximum native power grid G — F is a cactus and every edge of G — F that lies on 


11 



a cycle has infinite capacity. For any feasible flow / there exists an equivalent feasible flow 
f that is electrically feasible for Gp- 

Proof. We first construct an equivalent flow /' as follows. For each biconnected component 
C oi G — F that is a cycle, we consider the restriction fc of / to C. By Lemma |5.4[ there 
exists a unique flow ff, equivalent to fc that is electrically feasible for G. We now define 

w, , I f'riu.v) if u.v are in a cycle C 
]f{u,v) otherwise 


Note that changing in / the flow along the edges of a cycle G to the values determined by ffi 
preserves the net out-flow at every vertex, and hence f' is a flow equivalent to /. We claim 
that /' is electrically feasible. To see this, observe that each block of G — F is either a single 
edge or a cycle C. In the former case, /' is trivially feasible on the block. In the latter, we 
have that /' coincides on C with ff,, which is electrically feasible. By Corollary 5.2 /' is 
electrically feasible. □ 


Let ei,..., e/j be the edges of a cycle in Gp and fi be a flow on an edge in cycle 
C. We abbreviate the susceptance B{ei) on an edge in a cycle by Bi. The maximum 
susceptance is denoted by i?max with Bmax = maxi<i<fc(i?i) for all i = l,...,fc. The 
minimum susceptance Bmin is defined analogously. In practice, the requirement for infinite 
capacity in Theorem |5.5| is unnecessary. In fact, we can bound the sufficient large capacities 
of Theorem |5.5| by rearrange the equation of the proof of Lemma |5.4| such that the change 
of flow is bounded by the ratio of maximum to minimum susceptance times the average flow 
in the cycle G that is 

E fe 1 — D . U 

■ , ^min ^ 

1=1 Bi 

We refer back to Figure]^ which in addition to the previously mentioned parameters 
also shows the size of a minimum feedback vertex set with respect to cacti. In all cases the 
number of vertices for full control is between the size of feedback vertex sets with respect 
to forests and cacti. For the cases 14, 57 and 118, the minimum number of controllers for 
achieving full control indeed results in a native power grid that forms a cactus, although 
they do not necessarily achieve the smallest feedback number due to some influence of line 
capacities. 


6 Grid Operation Under Increasing Loads 

In the previous section we have seen that typically selecting a small fraction of the buses as 
flow control buses suffices to achieve full control in the network. In this section we study 
what happens when even fewer flow control buses are available and whether few flow control 
buses allow a better utilization of the existing infra-structure in the presence of increasing 
loads. 

To measure the controllability in the presence of very few flow control buses, we simulate 
a load increase by a factor p in the power grid by decreasing all line capacities by the factor 
1/p. This has the effect that the overall demand remains constant and thus any change 
of costs is due to flow redirections. It is then expected that, once the load increases, the 
network without flow control buses will require significantly higher operating costs, since the 
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Figure 4: Operation costs of case57 for OPF and the hybrid model with 1 and 2 control 
buses with respect to the load factor p. 


main criterion for determining the generator outputs becomes the overall feasibility of the 
flow rather than the cost-efficient generation of the energy. At some point, the load increases 
to a level where, by means of changing only the generator outputs, a feasible energy flow 
cannot be found. We compare the operation costs to solutions in power grids with a small 
number of flow control buses. Specifically, our plots show two things. First, the operation 
costs for various small numbers of flow control buses and, second, the operation costs and 
the number of flow control buses for achieving full control in the network with respect to 
the load increase factor p. 

Of course these operation costs again vary depending on the value of A. Since most 
related work ignores line losses, we consider only the case A = 1, i.e., only generation costs 
are taken into account. Varying A changes the objective value, but it does not influence the 
existence of solutions with a certain number of flow control buses. Recall from the plot in 
Figurej^that, if the load increase p is small, full control can be achieved without flow control 
buses for A = 1. In the IEEE instances all lines have very large capacities, often much larger 
than even the total demand in the network, e.g., 9900 MW in casel4 and case57. To better 
highlight the interesting parts, similar to the work by Lima et al. 15 , we first scale all line 


capacities such that the smallest capacity is equal to the total demand of the consumers as 
given in Table This changes neither the existence nor the cost of solutions. We increase 
the load until the flow model becomes infeasible; at this point a feasible solution cannot be 
achieved by adding flow control buses and adding additional lines to the network becomes 
unavoidable. 

Figure shows the results of our experiment for the power grid case57. To improve 
readability, all costs have been rescaled by the total demand in the network, and thus give 
the cost per MWh. The black curve shows the operation cost with sufficient control buses 
for full control. The dotted staircase curve shows the number of flow control buses that are 
necessary to achieve full control. Moreover, for each number of flow control buses from 1 up 
to the number required at the point when further load increase makes the instance infeasible, 
we show the optimal operation costs with this number of flow control buses. Finally, the 
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bold gray curve shows the operation cost with OPF, i.e., without any control buses. The 
plots for the other IEEE instances can be found in Appendix [C| 

As expected, increasing loads result in increasing operation costs. Interestingly, very few 
control buses suffice for extending the maximal feasible operation point. This is emphasized 
by the curve for two control buses in Figure which continues to a load increase of fac¬ 
tor 23.09, whereas OPF works only for up to an increase of roughly 17.27 and exhibits a 
significant increase in operation costs at higher loads. In contrast, when using flow control 
buses, the costs start increasing much later and more moderately. Interestingly, the solution 
with one control bus remains roughly equivalent to the solution with two control buses until 
shortly before the end of its feasibility range. This example shows that control buses indeed 
extend the feasible operation point and also decrease the corresponding operation costs even 
if there are only very few controllers available. 


7 Conclusion 

Assuming the existence of special buses that control the flow on all their incident trans¬ 
mission lines, we have presented a hybrid model for including some flow control buses. In 
this model, we have shown that relatively few control buses suffice for achieving full control. 
Further, we scaled the load of the network and showed that even fewer flow control buses 
improve the loadability and have a lower cost increase compared to OPF. 

Our work shows the benehts of augmenting power grids with flow control devices. Using 
our theoretical model, we were able to explain our empirical observations on controller 
placement with graph-theoretical means. While this also explains previous observations of 
Gerbex et al. , the main drawback is that the model is based on several strong, simplifying 
assumptions. 

Future work should consider more realistic power grid models both in terms of the control 
units, which are placed on transmission lines rather than buses, and using the AC power 
grid model. 
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A Problem Formulation 


In this appendix, we present the problem formulation from Sectionj^as an integer linear pro¬ 
gram (LP) formulation. In the LP below, we minimize the generation costs Cg and losses Ci 
shown in Equationunder flow and electrical constraints. The main flow constraints com¬ 
prise the conservation of flow, demand and generation constraints, and capacity constraints 
shown in Equations and respectively. Whereas the electrical constraints de¬ 

scribe the electrical feasibility shown in Equation Recall that each consumer u has a 
power demand ci„ € K and F C V is the set of flow control buses. 


min 

/ 


Ca(/) 

= >'-CgU) + {l-X)-ct{f) 


(11) 

s.t. 

E 

f{v,u) 

= 0 

u e E \ (Eg U Eg) 
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Figure 5: Trade-off of generator costs and costs of the losses depending as A varies from 0 
to 1. The square cross marks the solution computed by OFF. 
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Figure 6: Operation costs of case6 to casellS for OPF and the hybrid model with their 
control buses with respect to the load factor p. 
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